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ABSTRACT 

Quasi-coherent oscillations have been detected in the extreme ultraviolet flux 
of the dwarf nova SS Cygni during observations with the Extreme Ultraviolet 
Explorer satellite of the rise and plateau phases of an anomalous outburst in 
1993 August and a normal outburst in 1994 June/ July. On both occasions, the 
oscillation turned on during the rise to outburst and persisted throughout the 
observation. During the 1993 outburst, the period of the oscillation fell from 
9.3 s to 7.5 s over an interval of 4.4 days; during the 1994 outburst, the period 
fell from 8.9 s to 7.19 s (the shortest period ever observed in SS Cyg, or any 
other dwarf nova) within less than a day, and then rose to 8.0 s over an interval 
of 8.0 days. For both outbursts, the period P of the oscillation was observed 
to correlate with the 75-120 A count rate Jeuv according to P oc /guv 94 - ^ 
magnetospheric model is considered to reproduce this variation. It is found that 
an effective high-order multipole field is required, and that the field strength 
at the surface of the white dwarf is 0.1-1 MG. Such a field strength is at the 
lower extreme of those measured or inferred for bona fide magnetic cataclysmic 
variables. 

Subject headings: stars: individual (SS Cygni) — stars: magnetic fields — stars: 
novae, cataclysmic variables — stars: oscillations 



- 3- 



1. Introduction 

Rapid periodic oscillations are observed in the optical flux of high accretion rate 
("high-M") cataclysmic variables (CVs; nova-like variables and dwarf novae in outburst) 
( Patterson 19811 ; IWarner 19951 ; IWarner 19961) These oscillat ions have high coherence 



(Q ~ 10 4 -10 6 ), periods P ~ 10-30 s, amplitudes of less than 0.5%, and are sinusoidal to 
within the limits of measurement. They are referred to as "dwarf nova oscillations" (DNOs) 
to distinguish them from the apparently distinct longer period, low coherence (Q ~ 1-10) 
quasi-periodic oscillations (QPOs) of high-M CVs, and the longer period, high coherence 
(Q ~ 10 10 -10 12 ) oscillations of DQ Her stars. DNOs have never been detected in dwarf 
novae in quiescence, despite extensive searches; they appear on the rising branch of the 
dwarf nova outburst, typically persist through maximum, and disappear on the declining 
branch of the outburst. The period of the oscillation decreases on the rising branch and 
increases on the declining branch, but because the period reaches a minimum about one 
day after maximum optical flux, dwarf novae describe a loop in a plot of period versus 
optical flux. 

The dwarf nova SS Cygni routinely exhibits DNOs during outburst. Optical oscillations 
have been detected at various times with periods ranging from 8.2 s to 10.9 s ( Patterson 



ct al. 1975 ; Home fc Gomer 198U| ; Patterson 1981 ). During one outburst, the period was 



observed over a interval of ~ 6 days to fall from 7.5 s to 7.3 s and then rise to 8.5 s 
( |Hilde brand et al. 1981] ) . At soft X-ray energies (E rs 0.1-0.5 keV), oscillations have been 
detected in HE AO 1 LED 1 data at periods of ~ 9 s and 11 s ( |C6rdova et al. 1980| ; |C6rdova 



et al. 1984] ) and in EXOSAT LE data at periods between 7.4 s and 10.4 s QJones fc Watson 



1992| ). In this Letter, we describe observations with the Extreme Ultraviolet Explorer 
satellite (EUVE; |Bowyer fc Malina 199l |; IBowyer et al. 1994] ) of the EUV oscillations of 
SS Cyg. 
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2. Observations 

Target-of-opportunity observations of SS Cyg in outburst were made with EUVE in 
1993 August (MJD 9216.58 to 9223.12; MJD = JD - 2440000) and 1994 June/July (MJD 
9526.67 to 9529.78 and 9532.54 to 9536.94). On both occasions, SS Cyg was detected with 
the EUVE Deep Survey (DS) photometer and Short Wavelength (SW) spectrometer. The 
optical and DS count rate light curves of the outbursts are shown in Figure 1 of Mauche 
(1996a). On both occasions, the optical flux was above V = 10 for pa 16 days, but the 1993 
outburst was anomalous in that it took pa 5 days for the light curve to reach maximum, 
whereas typical outbursts (such as the 1994 outburst) reach maximum in 1 to 2 days. The 
EUV light curve of the 1993 outburst rose more quickly than the optical, and the rise of 
the EUV light curve of the 1994 outburst was delayed by ~ 1 day relative to the optical. 
On both occasions, the EUV spectrum evolved homologously over roughly two orders of 
magnitude in luminosity. The EUV spectrum is very complex, but it can be very crudely 
approximated by a blackbody with a temperature of kT pa 20-30 eV absorbed by a neutral 
hydrogen column density of Nr ~ 7-4 x 10 19 cm -2 ([Mauche et al. 1995Q . 

To detect oscillations in the EUV flux of SS Cyg, we proceeded as follows. First, we 
selected valid intervals when the source was above the Earth's limb and both the DS and 
SW instruments were in operation. Only those intervals longer than 600 s were retained for 
further consideration. For the 1993 (1994) observations, 89 (75) valid intervals were defined 
with a net time of 155 (103) kiloseconds. Second, we calculated the background-subtracted 
count rate in the DS and SW instruments for each of the valid time intervals. For the SW 
instrument, only those photons between 75 and 120 A were included in the sum; above 
and below these limits, the spectrum is dominated by noise QMauche et al. 1995|) . Third, 
for each of the valid time intervals, we constructed background-subtracted DS count rate 
light curves with 1 s time resolution. Fourth, we calculated the power spectra of these light 
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curves using the XRONOS software package. As with the HEAO-1 flCordova et al. 1980Q 
and EXOSAT ( |Jones fc Watson 1992|) soft X-ray power spectra, typical EUVE power 



spectra consist of a single line in a single frequency bin, with no power above the noise 
at any of the harmonics (Mauche 1996a| ). From this result, we understand that the EUV 
oscillation is (a) sinusoidal and (b) essentially coherent (Q = |AP/At| _1 > 6 x 10 4 ) over the 
~ 30 min valid time intervals. The period of the oscillation is then simply the inverse of the 
frequency of the bin of the power spectrum with the peak power, and its error is less than 
the bin width of P 2 /2048 ~ 0.03 s. Fifth, we calculated the amplitude of the oscillation 
during each valid time interval by phase-folding the data on the period appropriate to that 
interval and fitting a function of the form /(</>) = A + B sin(0 + <p ) where cf) = 2nt/P. The 
relative amplitude B/A was 14% ± 12% during the 1993 outburst and 16% ± 12% during 
the 1994 outburst, with a tendency in the 1994 data for the relative amplitude to decrease 
as the count rate increased. These amplitudes are comparable to those of the soft X-ray 
oscillations measured by HEAO-1 flCordova et al. 19801) and EXOSAT ([! ones & Watsonl 



The period of the EUV oscillation of SS Cyg detected by the above means is plotted as 
a function of time in Figure 1. For the 1993 outburst, the oscillation turned on during the 
rise to outburst on MJD 9218.77; its period was initially 9.31 s, fell over the next ~ 1.6 days 
to pa 7.6 s, and then asymptotically approached ~ 7.5 s over the next few days. For the 
1994 outburst, the oscillation turned on during the fast rise to outburst on MJD 9528.06; 
its period was initially 8.90 s, fell to 7.19 s (the shortest period ever observed in SS Cyg, or 
in any other dwarf nova) within less than a day, and rebounded to 7.42 s by the end of the 
observation on MJD 9529.77. When observations resumed 2.8 days later, the period of the 
oscillation was 7.59 s and rose slowly over the next few days to ~ 8.0 s. On both occasions, 
the oscillations switched on when the SW count rate was ~ 0.1 counts s _1 or when the DS 
count rate was ~ 1 count s -1 flMauche 1996a| ). For a coherent oscillation to be detected 
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below this count rate, its amplitude must be > 5%. 

Superposed on Figure 1 is the log of the 75-120 A SW count rate as a function of time. 
It is clear from this figure that the period of the oscillation anti-correlates with the SW 
count rate, being long when the count rate is low and short when the count rate is high. To 
quantify this trend, we plot in Figure 2 the log of the period of the oscillation as a function 
of the log of the SW count rate. The solid line connects the points chronologically starting 
in the upper left-hand corner. The upper panel of the figure shows the slow evolution of the 
period from 9.3 s to 7.5 s during the 1993 outburst. The lower panel shows the initial rapid 
evolution of the period during the 1994 outburst as the period fell from 8.9 s to 7.2 s and 
then rose to 8.0 s. Two aspects of this plot require immediate comment. 

First, despite the vastly different rate at which the period evolved in this diagram 
during the 1993 and 1994 outbursts, the slope and normalization of the trend of period 
with SW count rate is to first order the same. Using a function of the form P = Po^euv> 
where Ieuv is the 75-120 A SW count rate, the combined data is fit with parameters 
Po = 7.08 ± 0.03 s and a = 0.09385 ± 0.00003 if the errors on the period are set to produce 
xl = 1, but a more realistic estimate for the exponent is a = 0.094 ± 0.030. 

Second, there is no "hysteresis" in the trajectory of the 1994 outburst — it doubles back 
on itself instead of "looping," as do plots of the period of optical DNOs versus optical flux 
(see, e.g., Patterson 1981| ). This figure makes clear that to first order the period of the EUV 



DNOs of SS Cyg is determined solely by the SW count rate. Because the SW spectrum does 
not change during the observations, the bolometric correction for the EUV flux is fixed. 
We conclude that the period of the EUV DNOs of SS Cyg is determined by the EUV/soft 
X-ray luminosity, and, by inference, by the mass-accretion rate onto the white dwarf. 

Yet another conclusion follows from the above result that the period of the EUV 
DNOs of SS Cyg is a single-valued function of the EUV flux. Optical and EUV light curves 
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of SS Cyg ( Mauche 19961 ), U Gem ([Long et al. 1996|) , and VW Hyi QMauche 1996b| 
demonstrate not only the well-know result that the optical flux leads the EUV flux by ~ 1 
day on the rise to outburst, but that the optical flux lags the EUV flux on the decline from 
outburst. With the (as-yet unproven!) assumption that the period of the optical and EUV 
DNOs are equal, the looping trajectory of the period of optical DNOs versus optical flux 
is understood to be simply due to the time- dependent lead/lag of the optical flux relative 
to the EUV flux, and, by extension, the mass-accretion rate onto the white dwarf. [Quite 
happily, just such an effect is predicted by the accretion disk limit cycle mechanism for 
dwarf nova outbursts (see, e.g., Figure 13 of |Cannizzo 19931) .] This result, combined with 
the high coherence, amplitude, and luminosity of the EUV DNOs, strongly supports the 
long-held belief that the optical DNOs are produced by reprocessing of the EUV DNOs. 



3. Discussion 

Many models have been considered to explain the DNOs of high-M CVs ( |Warncr| 
|1995|) . The low period stability of the oscillations rules out the rotating white dwarf (the 
DQ Her mechanism) as well as non-radial pulsations of the white dwarf as the cause of the 
oscillations; pulsations are observed in high inclination systems, ruling out the eclipse by 
the white dwarf of luminous blobs of material in the inner disk and boundary layer; r-modes 
and trapped g-modes fail because more than one mode would be excited; oscillations of 
the accretion disk fail because they are not confined to a particular annulus and hence to 
a particular period. Viable mechanisms are more difficult to construct. Warner & Livio 
(1996) propose that the oscillations are due to the combined action of the differentially 
rotating surface layers of the white dwarf and magnetically controlled accretion. 

Despite any other indication that the white dwarf in SS Cyg is magnetic, we consider 
the requirements of a magnetospheric model for the DNOs of SS Cyg. Such a model 
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has been recently applied by Finger et al. (1996) to the quasi-periodic oscillations of the 
hard X-ray flux of the X-ray transient A0535+262. The Burst and Transient Source 
Experiment onboard the Compton Gamma-Ray Observatory observed a "giant" outburst of 
A0535 during 1994 February 3 through March 20. Finger et al. applied a magnetospheric 
model to the pulsar in A0535 to explain the observed variations in the QPO frequency 
i/qpo and spin-up rate of the neutron star. The rationale is a follows: the neutron star 
accretes through a disk; the disk will be disrupted by the magnetic field of the pulsar 
at some radius ro; the observed QPOs will be generated by some mechanism at either 
(a) the Keplerian frequency 2 , kvk{tq) = (GM*/?- 3 ,) 1 / 2 at r or (b) the beat frequency 
z/ bcat (r ) = ^x( r o) — u * between the Keplerian frequency at r and the spin frequency of the 
neutron star and torques on the magnetosphere interior and exterior to the corotation 
radius r co = [GM» / (27r^) 2 ] 1/73 will spin the neutron star up or down as ro varies with the 
mass-accretion rate. Assuming a dipole magnetic field B(r) = fi/r 3 where /x = B(R*)R% 
is the dipole moment, and determining ro by conservation of angular momentum for a 
Keplerian flow ( Ghosh fc Lamb 199 1| ), 

r = ir(GM,)- 1 V /7 M- 2/7 , (1) 

where K is a constant of order unity. Hence, vk^q) oc M 3 / 7 . With these relationships and a 
simple estimate for the accretion torque, Finger et al. could fit the observed variations in the 
QPO frequency and neutron star spin-up rate as a function of the measured 20-100 keV flux 
by identifying the QPO frequency with either vk{tq) or ^Wt^o)- Specifically, z/qpo oc 1^ 
over a variation of a factor of ~ 10 in the observed 20-100 keV hard X-ray flux Jhx- 
For subsequent reference, we note here that for the parameters derived by Finger et al., 
\i « 1 x 10 31 G cm 3 and, at the peak of the outburst, M ps 4 x 10~ 9 M© yr -1 . With these 
values, the above equation gives r ~ 1 X 10 9 cm wlx 10 3 neutron star radii. The strength 
of the magnetic field on the surface of the star is B(R*) kIx 10 13 G and its value at r is 
B(r ) w 1 x 10 4 G. 



- 9- 



Whereas z/q PO oc 1^ for A0535, z/ DNO oc Jew f° r SS Cyg. If a magnetospheric model 
applies to SS Cyg as it does to A0535, why is the variation of the DNO frequency such 
a weak function of the EUV flux? Variations in the bolometric correction relating the 
observed EUV flux to the EUV/soft X-ray luminosity and hence to M are ruled out because 
the EUV spectrum does not change throughout the outburst. Two alternatives exist, both 
of which appeal to different scalings of r with M than shown above: r oc M~ a where 
a ~ 0.29. The first alternative involves the accretion disk. Ghosh (1996) considered four 
standard disk models and found a = 0.15 for a one-temperature, optically thick, radiation 
pressure dominated disk; a = 0.25 for a one-temperature, optically thick, gas pressure 
dominated disk; and larger values of a for more exotic two-temperature optically thin disks. 
The data for A0535 are well fit by the scaling for the one-temperature, optically thick, gas 
pressure dominated disk: ro oc M _a25 , hence Vk{tq) oc M 0,38 . Although a one-temperature, 
optically thick, radiation pressure dominated disk gives a weaker scaling for Vk{tq) with M, 
it is very unlikely that such a model describes the disks of high-M CVs. 

If the structure of the inner disk of SS Cyg does not differ significantly from that of 
A0535, perhaps the structure of its magnetosphere does. The weak dependence of the 
frequency of the DNOs of SS Cyg on EUV flux suggests that the magnetosphere of SS Cyg 
is "stiffer" than the dipole magnetosphere of A0535. Assuming instead of a dipole magnetic 
field a single star-centered multipole field B(r) = rrii/r l+2 where m; = B^R^R 1 ^ 2 is the 
multipole moment ( [Arons 1993] ), 

r = K' (GM*y 1 ^ l+ Vmt /{4l+3) M- 2 ^ l+3 \ (2) 

hence, z/^-(r ) oc M 3 /( 4i + 3 ). To match the observed behavior of SS Cyg, we require I = 7l|. 
At the peak of the outburst, the EUV/ soft X-ray luminosity of SS Cyg is L « 2 x 10 33 erg s _1 
for kT = 20 eV QMauche et al. 1995Q , hence M = 4R*L/GM* w 3 x 10~ 10 M yr" 1 , 
where we have used M* = 1.2 M and i?* = 3.9 x 10 8 cm. Identifying i^dno 
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with z/^(r ), r = 5.9 x 10 8 cm pa 1.5 white dwarf radii (if instead we had used 
M* = 1.0 M , for which R* = 5.5 x 10 8 cm, r = 5.6 x 10 8 cm pa 1.0 white 
dwarf radii; the magnetosphere is crushed to the surface of the white dwarf). With 
these values and the above equation, the strength of the magnetic field at tq is 
B{tq) pa (GMj 1 / 4 M 1//2 r ~ 5//4 pa 5 x 10 3 G (independent of V.) and its value on the surface of 
the white dwarf is B{R*) = B(r ) (r /R*) l+2 pa 1 x 10 4 (1.5)' G ~ 2±{° x 10 5 G. 

A surface magnetic field strength of 0.1-1 MG is sufficiently low to be impossible to 
detect directly. In the presence of the light from the disk and secondary, field strengths below 
pa 3-5 MG are undetectable via circular polarization measurements in the optical or infrared 
( [Stockman et al. 1992 ). Direct detection of the cyclotron harmonics is difficult because the 



cyclotron fundamental is at A = 0.54(5/0.2 MG) mm. Zeeman splitting broadens spectral 
lines by a measurable AA/A = 1.2 x 10" 3 (A/6563 A) (5/0.2 MG), but Stark broadening 
is significant in the high-temperature, high-density photosphere of the white dwarf, and 
the disk contributes lines with Doppler widths AA/A = 8.4 x 10 _3 (r/10 10 cm) -1 / 2 sini. 
Furthermore, unlike some other dwarf novae (e.g., U Gem, VW Hyi), the white dwarf in 
SS Cyg is not observed above the disk and secondary during quiescence, even in the UV. 

Of more concern is the high order required of the multipole magnetic field. This comes 
about because of the weak dependence of the frequency of the DNOs of SS Cyg on EUV flux 
and, by inference, the slow variation of r with M. We have produced this slow variation 
by increasing the sensitivity of the magnetic field strength on radius: B(r) = mi/r l+2 with 
/ ~ 7 instead of the usual expression with 1 = 1. In general, the magnetic field can be 
expressed as an infinite sum of multipoles: B(r) = S^rnj/r' 4 " 2 = S / c ^ 1 5/(r). The behavior 
of SS Cyg requires that at tq the strength of the I ~ 7 multipole component of the magnetic 
field be significantly greater than the sum of all other components: B 7 (r Q ) ^> Bi(r ) for 
/ 7^ 7. This condition is quite restrictive, as it requires that the corresponding surface fields 
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Bi(R^) <C (i?*/r ) 7 ~' B 7 (R*) ~ (1.5)'~ 7 B 7 (R*). It is unlikely that such a strong inversion 
in the distribution of multipole components arises naturally — the topology of the currents 
required to produce such a field would be complex, to say the least. 

It is much more likely that the intrinsic magnetic field of the white dwarf of SS Cyg 
is a low-order multipole. The magnetic fields of AM Her stars are predominantly dipolar, 
though often the field strengths at the two poles are not equal, indicating that the dipole 
is off-center or equivalently that a few higher order terms enter into the expansion of the 
field. If the intrinsic magnetic field of SS Cyg is similarly predominantly dipolar, we need 
to explain why its magnetosphere behaves so much more "stiffly" than a dipole. A clue 
comes from the relative sizes of the magnetospheres of SS Cyg and A0535. In the latter, 
the magnetosphere extends to ps 1000 stellar radii at the peak of the outburst, while in the 
former it extends to only ps 1.5 stellar radii. Suppose that the accretion disk of SS Cyg is 
sufficiently diamagnetic as to exclude a significant fraction of the magnetic field of the white 
dwarf. During the outburst, the magnetosphere will be pinched inward by the accretion 
disk as the mass-accretion rate increases and the inner edge of the disk moves inward. 
Trapped in the orbital plane between the inner edge of the disk at r and the surface of the 
white dwarf at R*, the magnetosphere will become significantly distorted if the ratio r /-R* 
becomes small, as it does in SS Cyg. Expressing the distorted field as a sum of multipoles, 
it is clear that higher-order terms will become more and more important as the field is 
increasingly distorted. We do not imagine that the / ~ 7 term ever actually dominates the 
multipole expansion of the field, but that the magnetosphere comes to behave as if it did. 
In the context of the simply theory behind equation (1), we require simply that conditions 
conspire to produce B(r ) oc r ( -' +2 - ) with / ~ 7. Our derived value of B(r ) ps 5 x 10 3 G 
should still be a reasonable estimate, but the extrapolation of the strength of the field 
from r to the surface of the white dwarf is very uncertain. The simply scaling gives 
B(R*) ~ B(r ) (r /R*) l+2 « 1 x 10 4 (1.5)' G ~ 2±^° x 10 5 G for / = 7±|, but this can be 
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considered only an order-of-magnitude estimate. 



4. Conclusions 

If a magnetospheric model applies to SS Cyg, it appears likely that the surface magnetic 
field strength of its white dwarf is 0.1-1 MG. This range of values is orders of magnitude 
lower than the field strengths of AM Her stars (B m 10-80 MG; |Gropper 1990| ; [Beuermann 



fc Burwitz 1996), but may overlap with the field strengths of DQ Her stars (very uncertain, 



by B ~ 0.1-10 MG; Patterson 1994j) . With such a field strength, it is a challenge for SS Cyg 



in quiescence not to manifest the photometric variations associated with DQ Her stars; 
the limit on such variations in the optical is 0.001 mag (~ 0.01%) near 0.1 Hz (|Patterson| 



1995|) . In outburst, the accretion flow should be channeled down to the footpoints of the 
magnetic field and produce hard and soft X-rays in the manner of AM Her and DQ Her 
stars. However, neither the eponymous DQ Her nor V533 Her are hard X-ray sources 



( |CJ6rdova et al. 1981 ), demonstrating that magnetic accretion can take place without the 



production of hard X-rays. While SS Cyg in outburst is a known hard X-ray source ([Tones 



|fc Watson 1992| ; |Nousek et al. 199^ ; [Ponman et al. 1995[ ), it is not known to oscillate in 



hard X-rays; the limit is 6% when soft X-ray oscillations were observed by HE AO 1 ( [Swank 



1979|) . To produce a more stringent upper limit, XTE observations are required. Planned 
simultaneous observations of SS Cyg in outburst with XTE and EUVE will determine the 
shapes of the hard and soft X-ray light curves, the correlations between the soft and hard 
X-ray fluxes, and the extent to which the oscillations in the soft X-ray flux are manifest in 
the hard X-ray flux. 



The EUVE observations of SS Cyg could not have been accomplished without the 
efforts of the members, staff, and director, J. Mattei, of the American Association of 
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Variable Star Observers; EUVE Deputy Project Scientist Ron Oliversen; and the staffs of 
the Center for EUV Astrophysics (CEA), the EUVE Science Operations Center at CEA, 
and the Flight Operations Team at Goddard Space Flight Center. The author greatly 
benefited by conversations with P. Ghosh, J. Patterson, R. Robinson, and B. Warner. 
M. Finger, P. Ghosh, M. Livio, and B. Warner generously provided copies of material prior 
to publication. The manuscript benefited from the comments of an anonymous referee. 
This work was performed under the auspices of the U.S. Department of Energy by Lawrence 
Livermore National Laboratory under contract No. W-7405-Eng-48. 
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Fig. 1. — Log of the 75-120 A SW count rate (filled circles with error bars) and log of the 
period (open circles) as a function of time for the (a) 1993 August and (b) 1994 June/July 
outbursts of SS Cyg. 



Fig. 2. — Log of the period as a function of log of the 75-120 A SW count rate for the (a) 
1993 August and (b) 1994 June/July outbursts of SS Cyg. The dotted line is an unweighted 
fit to the combined data: P = 7.08/p T ° T v 94 s. 
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